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Abstract 

The nucleon-nucleon interaction arising from the exchange of three pions and the excitation 
of A(1232)-isobars in intermediate states is studied. Approximating the A-propagator by the 
inverse AN mass-splitting, analytical expressions are derived for the spectral-functions of the 
isoscalar and isovector central, spin-spin and tensor NN-potentials in momentum-space. A trans¬ 
lation of the spectral-functions into coordinate-space potentials reveals that the main effect of 
these specific exchange and excitation mechanisms is a repulsive isoscalar central NN-potential. 


1 Introduction and summary 

The interaction between two nucleons based on chiral effective held theory has been studied in great 
detail over the past two decades, for reviews see [Dial Em. Most calculations employ the effective 
chiral Lagrangian formulated in terms of pions and nucleons that are chirally coupled with each other 
and to external sources. The excitation of baryon and meson resonances is encoded in the low-energy 
constants of the chiral pion-nucleon interaction terms of higher order. Such a framework provides 
an accurate description of the empirical nucleon-nucleon phase shifts, if extended to sufficiently high 
order. At the present time, calculations have been carried out up to order N^LO in the chiral 
(low-momentum) expansion [5l [6] and the incorporation of dominant contributions at order N^LO 
is underway. Still, it can be argued the explicit inclusion the A(1232)-isobar, the most prominent 
resonance in nuclear physics, allows one to resum a certain class of important contributions and 
therefore may lead to an improved convergence. In such a phenomenological extension of baryon 
chiral perturbation theory the delta-nucleon mass splitting, A = 293 MeV, is counted as an additional 
small scale parameter, comparable to the typical momentum p or the pion mass m.,^. The hnite-range 
parts of the 27r-exchange NN-interaction as generated by the pertinent triangle and box diagrams 
with intermediate A(1232)-isobars have been calculated hrst in ref. [7] using the Feynman diagram 
technique. Furthermore, the Bochum-Bonn group mu has worked out the subleading isospin- 
conserving corrections as well as the isospin-breaking corrections provided by the chirally interacting 
vrAA-system. 

The purpose of the present paper is to analyze the two-nucleon interaction which arises from the 
exchange of three pions in combination with excitations of A(1232)-isobars in intermediate states. In 
this sense it represents an extension of earlier works on the 37r-exchange NN-interaction in refs. mm 
to chiral effective held theory with explicit A(1232)-isobar degrees of freedom. 

The present paper is organized as follows. In section 2, the techniques to the calculate spectral- 
functions (or imaginary parts) of two-loop 37r-exchange diagrams in the heavy nucleon mass limit are 
prepared. In order to facilitate an analytical treatment (of at least the angular part) of the involved 
37r-phase space integrals, the non-relativistic A-propagator will be approximated by the inverse 
mass-splitting A~^. A consequence of this approximation is that the pion-induced excitation and 
deexcitation of virtual A(1232)-isobars can be conveniently condensed into (symmetrized) 27t- and 
dvr-contact vertices with a nucleon. In section 3, analytical expressions are derived for the spectral- 
functions corresponding to the isoscalar and isovector central, spin-spin and tensor NN-potentials in 
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momentum-space. This is separately done for five different classes of 37r-exchange diagrams which are 
grouped together according to the number of intermediate A(1232)-isobar excitations. For each class 
the resulting NN-potentials in coordinate-space are displayed in the distance region 1 fm < r < 2 fm. 
As a summary one hnds that the main effect of these specihc exchange and excitation mechanisms is 
a repulsive isoscalar central NN-potential Vc{r). Moreover, the contributions to the isovector tensor 
potential Wt{t) tend to cancel each other, while other components come out very small anyway. 
This overall result is remarkable in view of the fact that previous calculations of the chiral three-pion 
exchange NN-interaction in refs. mm have lead to a vanishing isoscalar central potential. 


2 Preparation 

Let us start with recalling the form of the (static) nucleon-nucleon potential in momentum-space: 


Tnn = Vc{q) + Vs{q)ai-a2 + VT{q)ai-qa2-q 

+ {Wc{q) + Ws{q) ai-a 2 + IFr(g) ai-q a 2 -q}fi-f 2 , (1) 


where q = \q\ denotes the momentum transfer between the initial- and hnal-state nucleons. The 
subscripts C, S and T refer to the central, spin-spin and tensor-type components, each of which occurs 
in an isoscalar {Vc^s,t) and an isovector version iWc,s,T)- The sign-convention for T/vat is chosen such 
that the usual one-pion exchange gives: = (5'A/2/7r)^(?^^ + q‘^)~^- The occurring physical 

parameters are: the nucleon axial-vector coupling constant qa = 1.29, the pion decay constant 
/tt = 92.4 MeV, and the average pion mass = 138 MeV. 

We are interested in the finite-range parts of Taw that are generated by certain dvr-exchange dia¬ 
grams. For this purpose it is sufficient to calculate the imaginary parts ImVc, s,t(*/^) and ImlFc’,s,r(^/^), 
obtained by analytical continuation to time-like momentum transfer q = z/r, with /z > These 

imaginary parts serve also as the mass spectra entering a representation of the (local) coordinate- 
space potentials in the form of a continuous superposition of Yukawa functions: 
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with A the spectral cutoff, set to a maximal value of A = 1.5 GeV in ref. [5]. For the isovector 
potentials Wc,s,T{r) a completely analogous representation holds. Note that the tensor potentials 
Vr(^) and IFr(r) are accompanied by the usual tensor operator 3di-f (J 2 T — ai- 02 - 

Making use of the Cutkosky cutting rules, the imaginary parts entering eqs. (2,3,4) are calculated 
from the pertinent two-loop dvr-exchange diagrams as integrals of the iVY —)• dvr —)■ iVY transition 
amplitudes over the Lorentz-invariant three-pion phase-space. In the center-of-mass frame this four¬ 
dimensional phase-space integral includes an angular part of the form: 



dxdy 

a/1 — x‘^ — y'^ — z'^ + 2xyz 


( 5 ) 


with X = v-ki and y = v-k 2 the directional cosines of two pion-momenta ki^ 2 - The unit-vector 
V is introduced by the four-velocity u" = (0,zF) of the heavy nucleon in the considered t-channel 
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Figure 1: Contact-vertices of two and three pions with a nucleon arising from the excitation of 
intermediate A(1232)-isobars. 


kinematics NN^ Stt^NN HD]. The integration region E in eq.(5) is an ellipse x^+y‘^—2xyz < l—z‘^ 
with semi-axes \/l -f 2 ; and \/l — z, where 2 ; = /ci • ^2 (see also eq.(28)). In order to allow for 
an analytical treatment of the angular integral in eq.(5) for all two-loop Svr-exchange diagrams, we 
approximate the non-relativistic delta-propagator by the inverse delta-nucleon mass splitting A“^. In 
the case of a simple one-loop 27r-exchange triangle diagram this approximation leads to the inequality: 


— Ami 2 — Ami 

-:-> — arctan--- 

fiA fi 2A 


( 6 ) 


where the complete result stands on the right hand side. A comparison of the associated central 
potentials reveals that the approximation using A“^ leads to an overestimation by about 10 — 20%. 
Although the kinematical situation is more complex for Svr-exchange, one can expect that the ap¬ 
proximation using A“^ provides an upper bound for the NN-potentials with a similar error margin. 

Once the energy-dependence of the delta-propagator is neglected, one can combine the direct and 
crossed A(1232)-isobar excitation to a 27r-contact vertex as symbolized by a hlled square in the left 
diagram of Fig. 1. The corresponding transition matrix-element reads: 


— 


i9A 

4/.^A 


GabcTc ^-{qaX Qb) “ 45ab Qa ' Qb 


(7) 


where % and % denote both outgoing pion-momenta. In order to arrive at this form the relation 
TaT^ = {26ab — i^abcT'c)/^ for the isospin (and spin) transition operators and the coupling constant 
ratio g-wNA/g-wNN = 3/\/2 have been used. The strong vrAA-coupling constant follows from the 
Goldberger-Treiman relation as g-j^NN = gAM/f^^ = 13.1, with M = 939 MeV the nucleon mass. The 
delta-nucleon mass splitting A has the well-known value A = 293 MeV. 

Fig. 1 shows on the right also a diagram with a dvr-contact vertex which arises from an additional 
direct coupling of the pion to the A(1232)-isobar. Summing over the six permutations of the isospin- 
indices (a, b, c) the corresponding transition matrix-element reads: 


A^SttAA 


AOpA"^ I - 75eabc <fa- (<?6 X gj +qa-qb ^-qc (IS^aftTc - 76acrb - 7SbcTa) 

+qa-qc^-qb {786acTb - 76abTc - 76bcTa) + qb'qc^-qa {786bcTa - 76acTb - 76abTc) \ , (8) 


where qa, qb, qc denote outgoing pion-momenta. In order to £x the AAvr-vertex we use the coupling 
constant ratio g-wAA/g-wNN = 1/5 of the quark-model. Furthermore, the relation for the isospin 
(transition) operators, TaQbT^ = {7>i^abc — 6abTc + 4:6acTb — 6bcTa)/8, with 0;, the 4x4 A-isospin 
matrices, has been employed together with an analogous relation for the spin (transition) operators. 

While the decomposition into isoscalar and isovector contributions is obvious from the isospin 
factor V -|- fFri-f 2 of a 3vr-exchange diagram, a certain technique is needed in order to separate the 
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spin-spin and tensor-like components. First, one removes from a spin-dependent expression of the 
form ai-Aa 2 -B the Pauli-matrices cfi and a 2 - The spin-spin part ImV^ is then obtained by con¬ 
tracting the Lorentz-tensor A^Bj^ with the projector - 1 - and the combination 

— Ks) including the tensor part follows hnally by contracting A^Bp with the projector 
fi~‘^q°‘q^. After these manipulations the relevant expressions are given in terms of Lorentz-scalars, 
which can be easily translated from the s-channel {NN-^NN) into the f-channel {NN^NN). 


3 Calculation of spectral functions and r-space potentials 

In this section, we present the analytical expressions for the spectral-functions pertaining to the 
isoscalar and isovector central, spin-spin and tensor NN-potentials in momentum-space. This is done 
separately in hve subsections for the different classes of Svr-exchange diagrams, ordered according to 
the number of intermediate A(1232)-isobar excitations. In all cases the technical details are omitted 
and only the hnal results for the non-vanishing contributions are given. The calculated spectral- 
functions are then used to construct the corresponding NN-potentials in coordinate-space. These are 
presented in hve hgures displaying their dependence on the distance r in the region 1 fm < r < 2 fm. 


3.1 Single A-excitation of one nucleon 

We start with the Svr-exchange diagrams with a single A(1232)-excitation of one nucleon. The 
corresponding diagrams are compiled in ref. m, where they have been grouped into the classes X, 
XI, XII, XIII and XIV. The subleading 27r-contact vertex treated in ref. m includes M 2 nA as a special 
case, namely by setting the low-energy constants to the values: ci = 0, C 2 = —C 3 = 2 c 4 = g\/2A. 
Therefore, it suffices to evaluate the dtc-integrals in eqs.(7-22) of ref. [TT] for these parameter values. 
It is convenient to introduce the dimensionless variable u = /i/m,r > 3 and the auxiliary functions: 


R{u) = {u — l)\/{u — 3){u + 1), L{u) = In 


y/u — 3 + \/u + 1 


The contributions proportional to g\ which arise from the classes X, XI and XII read: 
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The additional contributions proportional to g\ from the classes XIII and XIV take the form: 
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Figure 2: NN-potentials arising from Svr-exchange with single A-excitation of one nucleon. 
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As indicated by the prefactor m^/A all the terms in eqs. (10-18) are of fourth power in small momenta. 
Thus they are counted of order N^LO in the chiral effective held theory with explicit A(1232)-isobars. 
Note also that the isoscalar central component ImVc vanishes identically due to the exact cancellation 
of the contributions from class XIII and XIV (see eq.(18) in ref. [TT]). 

The resulting local NN-potentials in coordinate-space are displayed by Fig. 2 in the distance 
region 1 fm < r < 2 fm. In order to obtain curves with less rapid decrease in r, we have divided out a 
Yukawa function exp(—3m7rr)/r with the decay-length (dm.^.)”^ = 0.48 fm. Such an unsealed Yukawa 
potential has a strength of (24.1, 5.65,1.49) MeV at distances r = (1.0,1.5, 2.0) fm, respectively. One 
observes from Fig. 2 that the strongest component is an attractive isovector tensor potential IFr(r), 
followed by a repulsive isovector central potential IFc'(^) of about half that magnitude. The spin-spin 
potentials, Vs{r) and Wsir), turn out to be particularly weak. 
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Figure 3: Diagrams related to Svr-exchange with single A-excitation of each nucleon. Upside-down 
reflected diagrams are not shown. Appropriate symmetry factors 1/2 are indicated. 

3.2 Single A-excitation of both nucleons 

Next, we come to the Svr-exchange diagrams with a single A(1232)-excitation of both nucleon. A 
representative set of two-loop diagrams with pion-bubbles and pions running in zigzag is shown 
in Fig. 3. The two-nucleon irreducible contributions coming from the planar and crossed bubble 
diagrams add up to the following results for the spectral-functions: 
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where A(w) = -f m^ -|- 1 — 2{w‘^v? + + m^) denotes the (kinematical) Kallen function. The 

dimensionless integration variable w introduced here is the invariant mass of a pion-pair, divided by 
the pion mass m^. 

The other diagrams in Fig. 3 with pions running in zigzag lead to contributions to the spectral- 
function which can be given in closed analytical form: 
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Figure 4: NN-potentials arising from 37r-exchange with single A-excitation of both nucleons. 
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The magnitudes of the pion-momenta /ci 2 , divided by and their scalar-product are given by: 
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where the last relation follows from energy conservation. The upper and lower limits of the U 2 - 
integration are = |(m —wii/civ^M^ — 2a;iM — 3/ \/u^ — 2uiu + 1), with ui lying in the interval 1 < 
oji < {v? — 3)/2m. Note that the non-algebraic result arcsin(x)/Vl — from the angular integral in 
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Figure 5: Diagrams related to 37r-exchange with double A-excitation of one nucleon. Left-right 
reflected diagrams are not shown. Appropriate symmetry factors, 1/6 and 1/2, are indicated. The 
third diagram including a Weinberg-Tomozawa 27r-contact vertex vanishes. 

eq.(5) prohibits a further step in the analytical integration over the 37r-phase space. In the chiral limit 

—)■ 0 the double-integrals in eqs.(25,26,27) have the values —M‘^7r^/2520, —M®(7r^/7 -|- 17/24)/360 
and m®( 31/24 — 7r^/7)/720, thus these spectral-functions become proportional to /r®/A^. 

Fig.4 shows the resulting NN-potentials in r-space, multiplied again with rexp(3m^r). One 
observes as the largest component a repulsive isovector tensor potential Wt{t), and for the hrst time, 
a repulsive isoscalar central potential Vc{t) of about half that strength. Note that the underlying 
spectral-functions, written in eqs.(22,23,24) have a rather simple analytical form. Other components 
of the NN-interaction turn out to be very small for this class of 37r-exchange diagrams. 


3.3 Double A-excitation of one nucleon 


Now we consider the diagrams in Fig. 5 including (on the left) a double A(1232)-excitation of one 
nucleon. The hrst two diagrams involving the leading order chiral dvr-contact vertex and the chiral 
vrvr-interaction give rise only to isovector spin-dependent contributions (proportional to g^) of the 
form: 
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Interestingly, the third diagram in Fig. 5 involving the Weinberg-Tomozawa 27r-contact vertex leads 
to vanishing spectral-functions \mWsgp = 0. This zero result is not obvious and it is obtained after 
the third step of the analytical integration over the pion-energy U 2 - Hence, one is left with the 
right diagram in Fig. 5 involving three ordinary vrA-couplings. The corresponding contributions are 
proportional to g\ and the double-integral representations of the spectral-functions read: 


ImVc = 


2(87r/2)3A2 


dujiduj 2 kik 2 \/1 — z'^ arccos(—x), 


(31) 


Z^<1 





























ImWs = 


(8^1^3A2 jj {^ + - 2ul -6)-‘^ + 2u,[u + U2{^ - u^)_ 

Z^<1 

zk r 1 ‘\ijo‘^ up' 

+ —^ uf{9u — AUJ 2 ) -^ H—^(9 + 47u 2 — 60u^ — 105u2u) + 5ui(2uj2 + 3m) 

6ki 5 5 


56 + 47a;2- 


56 + 47a;ol 9 , 9 9 „ ,1 ki, , 

-\-\bUJ 2 U --- + ^ {uJ2 — 1) +3) — U1i(u!2 + 3m)- 7r{zk2y 

5 J L5 J 3 

H- ^ ^ u^u{2u}2U — — 1) + uf(6u‘^ul — 4:U}2{2u^ + 3m) + m"^ 

kik2vl — z"^ I ^ 


+ 9m 5 j ^ ^ ^2a;2(4M^ + 15 m^ + 5) — m^ — 16m^ — 7u^ + 


u — 5 9 


lm{l/WT-Ws) = (3^J2^3A2 jj d^^duJ2 |^( 5 m -Ui- UJ 2){2 + UJ1UJ2) - 4 u^UJiUJ 2 

22<1 

+ —^ 5m(3m — UJ 2 ) + 5ui(Au2 — m) + a;?(12 — 15 m^ — 13m;2'W + 6a;9) 

3ki L 

+a;^(5M - 2^2) + 22:^a;i(a;2 - l)(i^i + 1^2 - 5m) - 

J OfC^ 

arccosf—2;) 9, , ^ ^ x , 1 

H- . = (1 -|- M ) (m — 2uji — 2CJ2) 4“ 4:UJ\UJ2U 

kik2vl — '- 

X ^(1 + M^) + cj^M + ^(6 c<; 2 M — 5 — 7m^) (33) 

We note as an aside that in the chiral limit m^r —)■ 0 the double-integrals in eqs. (31,32,33) have the 
values M^7r^/1260, M®(7r^ — 17/48)/540, and M®(7r^/7 -|- 31/480)/108. Hence, these spectral-functions 
become again proportional to p®/A^. It should be pointed out that the spin-isospin structure of 
the 37r-contact vertex Mst^aa in eq-(8) allows only for an isoscalar central component Vc and for 
isovector spin-spin and tensor components Ws,t- 

The resulting NN-potentials in r-space, multiplied with rexp(3m^r), are shown in Fig. 6. The 
dominant feature is now a repulsive isoscalar central potential Vd(r), whereas the other potentials 
hFs',r(^) are suppressed by about one order of magnitude. 

3.4 Double and single A-excitation of either one nucleon 

We continue with the 37r-exchange mechanisms accompanied by double and single A(1232)-excitation 
of either one nucleon. These are represented by the hrst two diagrams in Fig. 7 (and their mirror 
partners). In this case the pertinent spectral-functions can be given in closed analytical form and 
they read: 
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Figure 6: NN-potentials arising from 37r-exchange with double A-excitation of one nucleon. 
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Figure 7: The hrst two diagrams show Svr-exchange with double and single A-excitation of either 
one nucleon. The third diagram represents Svr-exchange with double A-excitation of both nucleons. 
The pertinent symmetry factors, 1/2 and 1/6, are indicated. 
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where the auxiliary functions R{u) and L{u) have been dehned in eq.(9). 

The resulting NN-potentials in r-space, multiplied with rexp(3m^r), are shown in Fig. 8. The 
dominant feature is again a repulsive isoscalar central potential Vc{r), while the other potentials 
are suppressed. 


3.5 Double A-excitation of both nucleons 


Finally, there is the Svr-exchange mechanism with double A(1232)-excitation of both nucleons. It 
is represented by the right diagram in Fig. 7. While the imaginary part of the isoscalar central 
NN-potential ImVd can be given in terms of the short expression: 
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Figure 8: NN-potential arising from 37r-exchange with single and double A-excitation of either one 
nucleon. 


rather lengthy formulas result from the multitude of spin-dependent terms in written in 

eq.(8). The imaginary parts of the isovector spin-spin and tensor NN-potentials \mWs,T read: 
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Figure 9: NN-potentials arising from 37r-exchange with double A-excitation of both nucleons. 

The corresponding NN-potentials in r-space, multiplied with rexp(3m^r), are shown in Fig. 9. 
The dominant feature is again a repulsive isoscalar central potential Vc{t). It is however smaller 
than those encountered in the previous two subsections. 

4 Conclusions 

In this work we have studied the two-nucleon interaction which arises from the exchange of three 
pions and the excitation of A(1232)-isobars in intermediate states. The pertinent two-loop spectral- 
functions have been calculated in the approximation of neglecting the energy-independence of the 
non-relativistic A(1232)-propagator. A subsequent analysis of the NN-potentials in r-space has 
revealed that the main effect of these sp^cihc exchange and excitation mechanisms is isoscalar central 
repulsion. Other components, such as IF'r(r) receive canceling contributions, or come out very small 
anyway. From the spectral-functions collected in section 3 one can reconstruct the NN-potentials in 
momentum-space via subtracted dispersion relations. This way one can include the interaction terms 
calculated in this work (in chiral effective held theory) into future hts to empirical NN phase shifts. 
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